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. Abstract 

A new algebraic Bethe ansatz scheme is proposed to diagonalise classes of in- 
tegrable models relevant to the description of Bose-Einstein condensation in dilute 
alkali gases. This is achieved by introducing the notion of Z-graded representations 
^ | of the Yang-Baxter algebra. 

1 Introduction 

(N 

Theoretical studies into the behaviour of Bose-Einstein condensates (BECs) continue at a 
prolific rate, motivated by the experimental successes of producing condensates of atomic 
alkali gases [J and superpositions of atomic-molecular alkali gases [|3], ^J. Many of the 
theoretical results to date have been obtained through use of the Gross-Pitaevskii mean- 
field theory and generalisations (e.g., see |5], |], |7|, |], |9|). However, such mean-field theory 
approaches have limited applicability in regions of the parameter space where quantum 
fluctuations dominate. In these cases, only an exact treatment of the model will give a 
reliable description of the physics. 

In our recent work we have shown that one model describing Josephson tunneling 
between two coupled BECs ]TIJ |TT[, and another that models coherent superpositions of 



atomic and molecular BECs fl2|| , can both in fact be solved exactly in the framework 
of the algebraic Bethe ansatz. Our intention here is to develop a new mathematical 
approach which allows us to extend this method to establish that very general classes of 
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Hamiltonians for BECs admit exact solutions. These classes of Hamiltonians cannot be 
solved in the usual form of the algebraic Bethe ansatz. 

In this Letter, three classes of solvable models relevant to Bose-Einstein condensates 
of dilute alkali gases are determined. This is achieved by formulating a new scheme for 
the algebraic Bethe ansatz by introducing the notion of Z-graded representations of the 
Yang-Baxter algebra. The first model we will present is a six-parameter generalisation of 
the canonical Josephson Hamiltonian (which can also be considered as a two site Bose- 
Hubbard model) describing a tunnel-coupled pair of trapped Bose-Einstein condensates. 
The Hamiltonian is also applicable to model solid state Josephson junctions and coupled 
Cooper pair boxes [13]. The second model we present has six free parameters and the 



third one has ten parameters. They both describe coherent coupling between atomic and 
diatomic molecular BEC's with additional interactions such as S'-wave scattering between 
the atoms, between the molecules, and between atoms and molecules. Such effects were 



not included in 12] but are important for a quantitative description of experiments 



rj|. Finally, we formulate the Slavnov formula for the scalar products between a Bethe 



eigenstate and an arbitrary Bethe vector, which facilitates the exact computation of form 
factors and correlations functions analogous to the results of [ITU . 



2 Bethe ansatz for Z-graded representations of the 
Yang-Baxter algebra 

The main ingredient in the study of exactly solvable quantum systems through the alge- 



braic Bethe ansatz [15, 16] is the Yang-Baxter equation 

Ru(u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u - v). 



Here Rjk(u) denotes the matrix in End (V <S> V <S> V) acting non-trivially on the j-th and 
k-th spaces and as the identity on the remaining space. The .R-matrix solution may be 
viewed as the structural constants for the Yang-Baxter algebra, denoted A, generated by 
the monodromy matrix T(u) 



R 12 (u - v)T x {u)T 2 {v) = T 2 (v)T 1 (u)R 12 {u - v). 



(2) 



The simplest case is that for the si (2) invariant i?-matrix, which will be the subject of 
our study, given by 
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with the rational functions b(u) = uj (u + rj) and c{u) = r]/(u + r/). 
Setting 

" A(u) B(u) 



(3) 



T(u) 



C(u) D(u) 



(4) 



2 



it follows from the defining relations (Q) that 



[A(u), A(v)} = [D(u), D(v)} = 0, 
[B(u), B(v)} = [C(u), C{v)\ = 0, 

A(u)C(v) = U ~ V + T] C{v)A{u) - -0—C(u)A{v), 
u — v u — V 

D(u)C(v) = U ~ V ~ V C(v)D(u) + -H—C(u)D(v). (5) 
u — v u — V 

Note that there are many more relations satisfied by the the generators of the Yang- 
Baxter algebra. However, those given above are the only ones needed for the algebraic 
Bethe ansatz procedure which we investigate below. For convenience, we extend A by a 
unit element / which we will represent by the identity matrix in any representation. 

We also introduce an auxiliary operator Z, called the grading operator, which satisfies 
the relations 

[Z,X(u)]=p{X(u)}.X(u), (6) 

where X = A, B, C or D and p{A(u)} = p{D(u)} = 0, p{B{u)} = 1 and p{C{u)} = -1. 
We call p{X{u)} 6 Z the gradation of X{u) and extend the gradation operation to the 
entire algebra by the requirement 

p {e.(f)} = p {6}+p{<f)} V0, <f>eA. 

This definition for the grading operator is consistent with the defining relations of A 
governed by (0). 

Let us now define a new class of representations of the Yang-Baxter algebra which 
we call Z-graded representations. We say that a vector space V, equipped with an endo- 
morphism z, is a Z-graded vector space, denoted (V, z), if it admits a decomposition into 
subspaces 

oo 

y= * 

k=—oo 

such that 

zV k = k.V k , k e Z. 

Note that some of the Vk may be trivial subspaces. Formally, the grading operator can 
be used to define the following projection operators 

J — — oo 

such that 

PkP\ = SklPh PkVj = SkjVk- 
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We say that a Z-graded vector space 

oo 

v= © K 

k=—oo 

is equivalent to V if for some j G Z there exist a vector space isomorphism between VJ. 
and Vj+k for all fc. This terminology is motivated by the fact that for a given (V, z) one can 
always generate another Z-graded space (V, z') through the mappings VI — > Vj+k, z' — > 
z — jl for any j G Z. 

For a given Z-graded V we say that it : A End V provides a Z-graded representation 
of A if 7r(Z) = z and the relations ar e preserved. In such a case we can write 

oo 

n(X(u))= XM) 

j=-oo 

where the matrices X(u,j) satisfy 

X(u,j)V k = iorj^k. 

More specifically, this means that for |^) e Vk we have 

tt(X(u)Y(v)) \i> k ) = X(u, k + p{Y(u)})Y(v, k) \ip h ) • 

In view of the equivalence of Z-graded vector spaces defined above, there can also exist 
equivalent representations. We can define a representation it' equivalent to tt by specifying 
some k G Z such that 

ir'{Z) = n(Z - kl) 

and for 



the matrices X'(u,j) are defined by 

X'(u,j) = X(u,j + k) Vj G Z. 

For any Z-graded representation it follows from ([5]) that the following hold: 

[A{u,j),A(v,j)\ = [D(u,j),D(v,j)} = 0, 
B(u,j)B(v,j-l) = B(v,j)B(u,j-l), 
C(u,j)C(v,j + l) = C(v,j)C(u,j + l), 

A(u,j)C(v,j + l) = U ~ V + r, C(v,j + l)A(u,j + l) 

u — V 

'' -C(u,j + l)A(v,j + l), 



u — V 



D(u,j)C(v,j + l) = '-C(v,j + l)D(u,j + l) 

u — V 



' -C(u,j + l)D(v,j + l), (8) 



u — V 
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From the defining relations (0) the transfer matrix defined by t(u) = A(u) + D(u) 
commutes for different values of the spectral parameter u; viz. 



[r(u),r(v)] = 0. 

Moreover, we may express the representation tt(t(u)) of the transfer matrix as 

oo 

n ( r ( u )) = Yl r ( u >ti 

j=-oo 

such that 

r(u,j)V k = Q for j ^ k 

and 

[t(uJ), r(v,k)] = Vj,k. 

Since p{r(u)} = 0, the diagonalisation of n(r(u)) is thus reduced to the diagonalisation 
of each of the matrices r(u,j) on the Z-graded component Vj, where we have 

[r(u,j), t(vJ)\ = 0. 

We may restrict our attention to the case of t(u, 0), as each r(u,j) is equivalent to some 
t'(u, 0) through the use of equivalent representations as introduced earlier. 

In order to formulate the algebraic Bethe ansatz solution for this class of representa- 
tions, we assume the existence of a pseudovacuum vector \x) € Vk such that 

A(u,k)\x) = a(u,k)\x) 

B{u,k)\ X ) = 

C(u,k)\ X ) + 

D(u,k)\ X ) = S(u,k)\ X ). 

The above implies that \x) is a maximal weight vector with respect to Z. Without loss of 
generality we can choose k = M, again due to the equivalence of representations discussed 
earlier, and look for Bethe states defined by 

= C(v u l)C(v 2 ,2)---C(v M ,M)\x). (9) 

It is easy to check that this Bethe state is symmetric with respect to the variables Vi, a 
feature which plays a crucial role below. Acting A(u, 0) and D(u,0) on the Bethe state 
we have 

M - + 

A(u, OMivi}) = a(u, M) H U y V *({ Vi }) 

i=l 
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M 



+ ^2Mi(u, • • • ,Vi-i,U,V i+1 , ■ ■ •,%), 

1=1 

M 

D(u, 0)tf -,%) = 5(u, M) TT U ~ Vi ~ V *({vi}) 

ti u ~ Vi 

M 

+ y^ j Af i (u, {vj})^(vi, ■ ■ ■ ,Vi-i,u,v i+1 , ■ ■ •,%), 



(10) 



with 



M - + 
= -^-a(^M)J] r ' + // 



U — Vj A Vj — V 
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M 



MM^}) = — — ^i,M)TT 

7/ — 11- J- J- 



U - V,; l>, - l>," 



(11) 



Requiring 

-Mi(u, {uj}) + {uj}) = 
forces \I/(i>i, • • - I'm) to be an eigenstate of t(u, 0) and leads to the Bethe ansatz equations 

«KM = TTW < = !,.. .,M. (12) 
The corresponding eigenvalue of the matrix r(u, 0) is 

M M 

A(«, W) = a(u, M) J] + ^ + S(u, M) J] M "^-^ . 

A U — Vi A - ± U — Vi 

i=l 4 i=l 1 

3 Explicit Z-graded realisations 

Next we give two nontrivial Z-graded realisations of the algebra A. One is expressible 
in terms of two Heisenberg algebras with generators dj, a\, i — 1, 2 and reads X(u,j) = 
X(u,j)P j with 

A(u, j) = u 2 + r]uN + rj 2 N 1 N 2 - r](N 1 - N 2 )uo(N + j7) 

-u 2 {N + ] I) + a\a l) 
B(uJ) = (u + u(N + jI) +r ] N 1 )a 2 + ^ 1 a 1 , 

£>(u,j) = a\a 2 + V - 2 . (13) 
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Above, Pj are the projections defined by ([?[), Ni = aja^, N = Ni + N 2 and uj(x) is an 
arbitrary polynomial function of x. Note that in the case when u(x) is constant, the 
above realisation reduces to that discussed in jn| [Tl[] and is factorizable into two local 
representations of the Yang-Baxter algebra expressible in terms of the two Heisenberg 
algebras. It is important to note that for generic u(x) no such factorisation exists. 

The representation acts on the infinite dimensional Fock space spanned by the vectors 

\m,n) = (a{) m (4) n |0> , m, n = 0, 1, 2, oo. (14) 

For this representation, we choose the pseudovacuum \x) as the Fock vacuum |0). The 
representation of the grading operator Z is chosen to be 

tt(Z) = M.I - N. 

We then have 

a(u,M) = u 2 -lu 2 (M), 5(u,M)=r]- 2 (15) 
and the Bethe ansatz equations become 

M 

^ 2 -^(M)) = l p~^T ( 16 ) 

for the diagonalisation of the matrix t(u, 0). The eigenstates (||) in this instance are also 
eigenstates of the total particle number N with eigenvalue M. 

Another Z-graded realisation of the Yang-Baxter algebra is X(u,j) = X(u,j)Pj with 

A(u,j) = -7]u 2 + u(l-7] 2 (K z + N c )-r]u(K z + N c + jI)) 

+r)K z - V 2 K z lu(K z + N C + jl) - V 3 N C K Z + V 2 cK + , 
B(u,j) = r](l-r]u-r]uj(K z + N c + jI)-r] 2 N c )K^-r]c(u-r]K z ), 
C(u,j) = r)c\u + T]K z ) -r)K + , 

D(u,j) = u-r]K z + rfc ] K_. (17) 

Above, the operators c, form a Heisenberg algebra, with iV c = c^c, and the operators 
K z , K + , K_ satisfy the relations of the su(l, 1) algebra [K z , K±] = ±K±, [K + , K_] = 
—2K Z . As in the previous example, u(x) is an arbitrary polynomial function of x and the 
above realisation is factorisable only in the case when u(x) is constant. 

For this representation, we choose the pseudovacuum \x) as the tensor product of the 
Fock vacuum |0) with a lowest weight state for the algebra su(l, 1) of weight k. The 
representation of the grading operator may be chosen as 

n(Z) = M.I -K t - N c . 

Then, 

a(u,M) = (1 -r)u-r)uj(M))(u + r)K), 6(u,M)=u-r)K (18) 
and the Bethe ansatz equations are 

(1- jjVi -r]u(M))[ =11 -—■ !9 
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4 Three models of Bose-Einstein condensates 



4.1 Model 1: Two coupled BECs 

Consider the following general Hamiltonian describing Josephson tunneling between two 
coupled Bose-Einstein condensates 

H = U u Nf + U 12 N 1 N 2 + U 22 Nl + n l N 1 + ^ 2 N 2 

_^( a t G2 + a t ai ). (20) 

The above Hamiltonian generalises the canonical Josephson Hamiltonian studied in [ITD, 



11] in that the couplings Un, U 22 for the S-wave scattering terms can be chosen arbitrarily. 



It also describes a pair of Cooper pair boxes with capacitive coupling [0. In the limit 
U22 — > 0, then (N 2 ) » (Ni), which can be considered as a single Cooper pair box coupled 
to a reservoir. 

It is an algebraic exercise to show that the Hamiltonian is related with the matrix 
f 0, 0) = A(u, 0) + D(u, 0) through 

H = [f(0,0) -r]- 2 + (aN + f3) 2 -r]aN -7]5N 2 ] . 

Here we have chosen u(N) = aN + (3 and the coupling constants are identified as 

2 2{U n + U22 - U 12 ) 



V 



a 



P 



a 



Noting that 



U11 - U22 

Hi + 1^2 

Uu + U22 



iV = r/" 1 ^(0,0), 
du 



the above demonstrates that the Hamiltonian ( p0|) is expressible solely in terms of the 
matrix f(u, 0) and its derivative. 

Since [H, N] = 0, the Hamiltonian is block diagonal on the Fock basis QI%P- Thus 
on a subspace of the Fock space with fixed particle number N, the diagonalisation of 
f(u, 0) is equivalent to the diagonalisation of t(u, 0) presented earlier in the Bethe ansatz 
framework. We then deduce that the solution of (BO) for the energy spectrum is 



e = Jl 

2 



N N 

- 2 n^-(^ + /3) 2 n— 

±J : Vi . 4" Vi 

i=i i=i 
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-7]- 2 + (aN + f3f - 7]aN - 7]SN 2 ] (21) 
where the parameters {v «} are subject to the Bethe ansatz equations 



N 

Vi -Vj-T] 



, . Vi - Vj + 7] 

4.2 Model 2: Homo- atomic- molecular BECs 

Next we turn our attention to a two-mode model for an atomic-molecular Bose-Einstein 
condensate with identical atoms. The Hamiltonian takes the form 

H = U aa N 2 a +U ac N a N c + U cc N 2 + fi a N a + fi c N c 

+VL(a ] a) c + c ] aa) (22) 

which acts on a basis of Fock states analogous to (|14]). For the case of 87 Rb, all of these 
parameters have been estimated from experiment (see ||). In the experiment described 
in [Q, the parameter U aa was varied significantly with a magnetic field. 

The Hamiltonian commutes with the total atom number iV = iV a + 2N C . In terms of 
a realisation of the algebra su(l, 1) through 

K+ = W t K _ = * Kz = w (23) 

+ 2 2 4 v ; 

one may establish the relation between the Hamiltonian and the corresponding transfer 
matrix f(u, 0) = A(u, 0) + D(u, 0) arising from the the realisation ( ]T7| ) of the Yang-Baxter 
algebra is 

H = a + 5(N/2 + 1/4) + 7 (iV/2 + 1/4) 2 + 2^ 2 fif (0, 0), 

with 



(if 

— (0, 0) = 2 - 77(77 + a)(iV/2 + 1/4) - rj(3. 
du 



Above we have chosen 



u{K z + N c ) = a (K z + N c )+P 
= a(N/2 + l/4)+(3 

and the following identification has been made for the coupling constants 

4U aa + U cc - 2U ac 

V = 2h ' 

U C c 4C/ aa 

a = — , 

2Q 

2/l c ~ 4/J a + AUaa ~ Vac 
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u„ 



a 



S 

7 



2/i c - £/a 



By the same argument as before, we conclude that the exact solution for the energy 
spectrum of (p2|) is determined by 

E = o + 5(M + k) +7(M + k) 2 



+2ti~ 1 kQ 



M M 

Vi — Tj T-r Vi + 7] 



(l-r)(a(M + K)+(3))H^—L-ll 



i=l 



(24) 



where the parameters Vi satisfy the Bethe ansatz equations 



[1 -rp)i -r)(a(M + K) + 



Vi + T)K 
Vi - T)K 



M 

n 



Vi -vj-rj 
Vi-Vj + rj' 



(25) 



For the representation ( P3[ ) of the su(l, 1) algebra there are two lowest weight vectors; 
viz. the Fock vacuum |0) and the one particle state (V |0). It follows that the allowed 
values for k in (^,^) are k = 1/4, 3/4. This demonstrates that the solution of the model 
depends on whether the total particle number iV = 2M + 2k — 1/2 is even or odd, the 
effects of which on the energy spectrum can be seen through numerical analysis (cf . |12j ) . 



4.3 Model 3: Hetero-atomic-molecular BECs 

The previous construction can be extended to model an atomic-molecular Bose-Einstein 
condensate with two distinct species of atoms, denoted a and h. For this case the Hamil- 
tonian takes the form 

H = U aa N 2 a +U bb Nt + U cc N 2 c +U ab N a N b + U ac N a N c + U bc N b N c 

+VaN a + fi b N b + fi c N c + VL(a)b ] c + c ] ba) (26) 

which commutes with the total atom number N = N a + N b + 2N C and I = N a — N b . Here 
the model acts on the Fock space spanned by the vectors 

|/ ) m,n> = (oWHct)»|o). 

In order to show the solvability of this model, we adopt the realisation of the su(l, 1) 
algebra given by 

K + = a)b\ K_ = ab } K z = ^H±3i±l ) (27) 

and observe that the operator X commutes with the su(l, 1) algebra in this representation, 
hence taking a constant value in any irreducible representation. Due to the symmetry 
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upon interchanging the labels a and b, we can assume without loss of generality that the 
eigenvalues of X are non-negative. In particular, note then that the lowest weight states 
for this realisation are of the form 

,t\m 



m ) = ( a T) m |0), m = 0,l,2, 



and K z \m) = (ra/2 + 1/2) \m). We conclude that the lowest weight labels k can be taken 
from the set {1/2, 1, 3/2,....} and the eigenvalue of I on the irreducible representation 
labelled by k is 2k — 1. 

For this case the relation between the Hamiltonian and the corresponding matrix 
f(u, 0) is 

H = a + 8{N/2 + 1/2) + A(iV/2 + 1/2) 2 

+ pi + vZ 2 + £J(iV/2 + 1/2) + r/~ 2 fif (0, 0) (28) 

with 

df 

— (0, 0) = 2 - 77(77 + a)(N/2 + 1/2) - V f3I - 777. 
Above we have chosen 

u(K z + N b ) = a(K z + N c ) + (3{2k — 1) + 7 
= a(iV/2 + l/2)+/3J + 7 

and the coupling constants are related through the relations 

Uaa + Ubb + U cc + U a b — U ac — Ubc 



77 = 



a 



Q 

f^cc — C4a — Ubb — U a b 

2U bb - 2U aa + U ac - U bc 
2Q ' 
2U aa + 2U bb + 2U ab - Ugc ~ Ubc + 2pc ~ 2/i a - 2jx h 

2tt 

Uaa + Ubb + Uab ~ 2/X a ~ 2fl b 

4 

2fJ, c — U ac — Ubc 



(3 

7 
a 

s 2 ■ 

A — Ucc, 

Ubb — U aa + Pa — Pb 

P = 2 ' 

U aa + Ubb — U a b 



V 



4 

U ac — U hc 
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The exact solution in this instance reads 
E 



o + 6{M + k) + \{M + k) 2 
+ p(2n - 1) + i/(2k - l) 2 + £(2k - 1)(M + re) 



M 



1 - r/(a(M + k) + /3(2k - 1) + 7)) Yl 



Vi-T] 



i=l 



M 

n 



n 



?Jj -Vj-T) 
Vi -Vj + T) 



(29) 



where the parameters Uj satisfy the Bethe ansatz equations 
[1 -77^-77 (a(M + k) + (3{2k - 1) + 7)] (m±HH 
and the total atom number is given by A^ = 2M + 2k — 1. 

5 Wave function scalar products 

Recall that in the usual algebraic Bethe ansatz for the algebra A there is a formula 
originally due to Slavnov |L7| (see also |L5|, |18|1 ) for the wave function scalar products. 
The Slavnov formula still applies in the Z-graded case and takes the usual form 



SM({Uj},{v k }) 



$({ % })*(H}) 
= *({«*})¥({«,■}) 
detT({itj},{t; fc }) 
detV({ Uj }, {v k }) 

with the entries of the M x M matrices T and V given by 



(30) 



T, 



ah 



d_ 

dv a 



A(it ft , {v k }),V a 



ab 



Ub - V a 



-,a,b = 1, 



M 



$({wj}) is the left vector 

$(m, • ■ ■ , u M ) = (x\ B(u M , M) ■ ■ ■ B(u u 1). 

Above, we have adopted the usual convention to scale the Yang-Baxter algebra such that 
5(u, M) = 1. Also, {ffc} provide a solution of the Bethe ansatz equation ( P^) and the 
parameters {uj} can be chosen arbitrarily. 

The Yang-Baxter algebra A admits a conjugation operation f : A — > A defined by 



A( U y = A(u), B(uY = C(u), C( U y = B(u), D( U y = D\ 



u 
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and extended to all of A through 



(0.0)t = 0t.0t ? v0, <peA 

such that the defining relations @ are preserved. Consequently the right vector $(t>i, • • • , vm) 
is also an eigenvector of the transfer matrix whenever the Bethe ansatz equations for the 
parameters {vt} are satisfied. However, it is apparent that the Z-graded representations 
( |I3| , [r7D we have introduced are not unitary, and generally 

$(K})t ^ 

On the other hand, numerical analysis we have undertaken for the above models indicates 
that for fixed particle numbers, and generic values of the coupling parameters, the energy 
spectrum is free of degeneracies. This is presumably due to the fact that the only Lie 
algebra symmetries for these models are u(l) invariances corresponding to conservation 
of particle numbers, and the non-degenerate spectra are examples of Hund's non-crossing 
rule fO, 1201. Whenever this is the case, we can conclude that 



$({^}) t = KV({vi}) 

for some constant K and the Slavnov formula can still be invoked for the computation 



of form factors and correlation functions (cf. the example of [10| where it was found 
K = ±l). 



6 Conclusion 

In conclusion we have introduced a new scheme for the algebraic Bethe ansatz to diag- 
onalise three classes of integrable models relevant to Bose-Einstein condensates of dilute 
alkali gases. The extension of this construction to other types of models, such as the 
Jaynes-Cummings model | 2T|| , and generalised Tavis-Cummings model discussed in |22| . 
is straightforward. 



This work is supported by the Australian Research Council. 
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